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ABSTRACT 

The Type lib SN 2011dh exploded in the nearby galaxy M51 (the Whirlpool Galaxy) 
and provides us with one of the best laboratory to study early high energy emission 
from SNe. We give here a comprehensive view of the X-ray properties of SN 201 ldh 
from the analyses of two pointed XMM-Newton early observations as well as of the 
full Swift X-ray Telescope (XRT) datasct (163 ks). Due to the high XMM-Newton 
throughput, we were able to satisfactorily fit the X-ray spectrum with two hot diffuse 
gas components including an additional absorption component to our Galaxy. A power 
law model provided a worse description of the data. In addition, the early Swift XRT 
light curve hints of a flux excess at early times ( £ 3 d) , consistent with the adiabatic 
cooling of stellar's photosphere a few days after the shock breakout. 
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1 INTRODUCTION 

Massive stars (M £ 8 M Q ) end their lives with a (core- 
collapse) supernova (SN) explosion. It is notoriously diffi- 
cult to detect X-ray emission from SNe. Dominant X-ray 
emission mechanisms comprise an energetic shock break out 
soon after the explosion (Campana et al. 2006; Soderberg et 
al. 2008) or the interaction of the shock front arising from 
the SN explosion with the circumstellar medium (CSM). In 
particular, collisions of the ejected material with the CSM 
(made by material ejected from the progenitor prior to ex- 
plosion) can heat up matter to T ~ 10 9 K, leading to X-ray 
emission. A reverse shock will also form propagating inward 
with a temperature T ~ 10 7 K. Alternatively, X-ray emis- 
sion can arise from Compton down-scattered gamma-rays 
of radioactive decay in the ejecta, or from inverse Comp- 
ton scattering of relativistic electrons on photospheric UV 
photons produced during the outburst. 

Nearby SNe represent an ideal target for X-ray studies. 
On May 31.893, 2011 UT a SN was discovered by several 
groups in the Whirlpool galaxy M51 (d ~ 8.4 Mpc, Feld- 
meier, Ciardullo & Jacoby 1997). The SN was classified as 
a Type lib (Arcavi et al. 2011; Marion et al. 2011) and 
named SN 201 ldh (Griga et al. 2011; Silverman, Filippenko 
& Cenko 2011). Given its closeness attempts to identify a 
progenitor in pre-explosion Hubble Space Telescope (HST) 
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images have been carried out. A possible progenitor has been 
identified even if it is not certain (Maund et al. 2011; van 
Dyk et al. 2011). Photometry of this object is consistent 
with a yellow super-giant star. Along these lines Bersten et 
al. (2012) modeled the optical light curve and the photo- 
spheric velocity evolution based on a set of hydrodynamical 
models, finding a very large progenitor radius of R„ ~ 10 13 
cm. A compact progenitor was instead suggested by other 
authors (Soderberg et al. 2012; Arcavi et al. 2011). In this 
case the identification of the progenitor was questioned and 
the radius estimate was based on radio properties leading to 
a radius estimate of R„ ~ 10 11 cm. 

SNe lib may be divided into two sub-classes based on 
the radius of the progenitor and on its mass loss history 
(Chevalier & Soderberg 2010). Compact progenitors give 
rise to fast shock waves (~ 0.1 c) and modulated radio light 
curves (SN cllb); extended progenitors (R„ ~ 10 13 cm) 
are characterized by slower shock waves and smoother ra- 
dio curves (SN ellb). Hydrodynamical light curve modeling 
gives strong evidence that SN 201 ldh belongs to the class 
of extended progenitors (R ~ 200 R@; Bersten et al. 2012). 

Here we report on X-ray observations on SN 201 ldh. 
We will first report on two XMM-Newton observations 7 
and 11 d after the SN explosion (Section 2). In light of these 
observations we will reanalyse Swift X-ray Telescope (XRT) 
observations (Section 3), expanding the results reported in 
Soderberg et al. (2012). Section 4 is dedicated to discussion 
and conclusions. 
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Table 1. XMM-Ncwton observations log 



Obs.l 


Exp. time 


Filt. exp. time 


Collected 


Obs.2 


Exp. time 


Filt. exp. time 


Collected 


Jun 7.2 2011 


(ks) 


(ks) 


photons* 


Jun 11.2 2011 


(ks) 


(ks) 


photons* 


pn 


9.8 


9.8 


719 (75%) 


pn 


9.7 


2.6 


123 (67%) 


MOS1 


10.6 


10.6 


165 (79%) 


MOS1 


6.4 


5.5 


80 (75%) 


MOS2 


10.7 


10.7 


181 (67%) 


MOS2 


6.3 


5.5 


60 (76%) 



* Photons within the extraction region in the 0.2—12 keV energy range for the pn and 0.3-10 keV for the MOS. In parenthesis we also 
reported the contribution of the source with respect to the estimated background. 




Figure 1. Left: smoothed image of all the Swift XRT observations (163 ks). Right ascension on the x-axis and declination on the y-axis. 
Right: smoothed image of the first XMM-Newton pn observation (~ 10 ks). Images were smoothed with a Gaussian function with 3 pixel 
kernel radius. SN2011dh is indicated in both images with a circle, that is also the extraction region adopted. 



2 XMM-NEWTON OBSERVATIONS 

XMM-Newton observed M51 twice based on a DDT program 
(P.I. Campana). The first observation took place on June 
7.207, 2011 UT and lasted ~ 10 ks (see Table 1). The thin 
filter was used on all the X-ray instruments for all the ob- 
servations. During the entire observing time the background 
remained low and we retained all the observing time. The 
second observation took place on June 11.196, 2011 UT for 
~ 10 ks (see Table 1). In this case a strong background flare 
affected the observation. We filtered the MOS observations 
requiring less than 10 c s _1 on the entire instrument. For 
the pn we used a 30 c s -1 threshold. The resulting exposure 
times are reported in Table 1. 

Data from the first observation were extracted from a 
500 pixel (25 arcsec) radius centered on source. The back- 
ground was evaluated in two nearby regions both with 800 
pixel radius on the same CCD surrounding the SN 2011dh 
location and free of sources and bad columns (see Fig. [T]). 
Chandra images show that within this extraction region 
there is just a faint source to the south-west (Soderberg et 
al. 2012). Given the brightness of SN 201 ldh, contamination 
of the data is minimal and does not affect our analysis. We 
collected 1065 and 263 photons for the first and second ob- 



servation, respectively. Data were group to have at least 20 
counts per energy bin and rmf and arf files were generated 
using the latest calibration products. 

We fit together the six spectra, keeping the same ab- 
sorption pattern for all of them (two TBABS models, one fixed 
to the Galactic value of 1.8x 10 20 cm -2 ; Kalberla et al. 2005) 
and allowing for the same normalization constant for MOS1 
and MOS2 in the two observations, while keeping the pn nor- 
malization fixed to one. The second spectrum is softer than 
the first one. We tried a few single component models. The 
null hypothesis probabilities are never larger than a few per- 
cent, indicating poor fits. We also tried a two more complex 
models. A model consisting of a power law and an emission 
from hot diffuse gas (MEKAL) provides an improvement (3% 
according to the F-test) but the resulting intrinsic column 
density is extremely large and not consistent with optical 
data (Arcavi et al. 2011; Ritchey & Wallerstein 2012, see 
also below). A model with two hot thermal gas components 
gives a better fit to the data (null hypothesis probability of 
15%) and an improvement over a single MEKAL of 2.7 a (0.7%, 
according to the F-test, Fig. [2J . 

Given the two best fit models we tried to fit the data 
without the presence of an intrinsic absorption component 
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Figure 2. XMM-Newton spectra of SN2011dh. Squares refer to 
pn data, circles to MOS1 data, and triangles to MOS2 data. Filled 
markers indicate data from the first observation, open markers 
from the second observation (black: 1 st pn, red: 1 st MOS1, green: 



1 st MOS2; blue: 2 n 
MOS2). 



pn, light blue: 2 nd MOS1, magenta: 2 n 



and evaluate the improve in the fit due to its inclusion with 
an F-test. In both cases the probability of a chance improve- 
ment is 2% (2.3 a). We also note that we assumed for the 
Galactic column density contribution a value of 1.8 x 10 20 
cm" 2 (Kalberla et al. 2005) which is larger than the val- 
ues derived from the Dickey & Lockman (1990) survey of 
1.6 x 10 20 cm" 2 or from the Stark et al. (1992) survey of 
1.5 x 10 20 cm -2 . This indicate that the estimated signifi- 
cance of our intrinsic component is a tight lower limit. 

In all models a non-negligible absorption component in 
addition to the Galactic one is needed. This indicates the 
existence a small but non-zero intrinsic absorption, possibly 
arising from the absorbing shell between the forward and 
reverse shock. The number of collected photons is too low 
to assess if a variation of the column density is occurring 
among the two observations. 

We adopted this model as our reference model (last row 
in Table 2). In this case the unabsorbed 0.5-10 keV flux is 
2.4 x 10" 13 erg cm -2 s" 1 and 1.1 x 10~ 13 erg cm -2 s -1 , 
respectively. These correspond to a 0.5-10 keV luminosities 
of 2.0 x 10 39 erg s" 1 and 9.3 x 10 38 erg s~\ respectively 
(adopting a power law, fluxes and luminosities are ~ 10% 
lower) . 



3 SWIFT X-RAY TELESCOPE 
OBSERVATIONS 

Swift promptly detected SN 2011dh about 3 days after the 
explosion (Margutti & Soderberg 2011). Swift observed M51 
112 times. 40 observations were carried out before the SN 
201 ldh explosion. Soderberg et al. (2012) provided a 3 a up- 
per limit at the SN position of 5.6 x 10 -4 c s" 1 (converting to 
a 0.5-10 keV luminosity of ~ 2 x 10 38 erg s" 1 ). The remain- 
ing 72 observations were targeted on SN 2011dh spanning 
the time period June 3 - August 15, 2011 for a total observ- 
ing time of 163 ks (Fig. [I]). The light curve was built using 



the UKSSDC light curve generatoir paying particular care 
not to include nearby sources (Fig. 36. Given the variable 
XRT background caused by CCD temperature variations, 
we prefer not to rely on pre-explosions observations. Each 
point in the light curve has a significance of 6 a. We assumed 
an explosion data of May 31.6, 2011 UT, based on the analy- 
sis of Soderberg et al. (2012). This sets the origin of the time 
axis in the Fig. 3. The overall curve can be fit with a power 
law decay with index —0.59 ± 0.16 (including the small un- 
certainty involved by the determination of the SN explosion 
of 0.2 d, see Fig. 1). The fit is acceptable with a reduced 
X 2 c d = 1-2 (20 d.o.f.). One can notice a flux excess at the 
beginning of the observing time, possibly related a SN shock 
break out (Campana et al. 2006; Soderberg et al. 2008). At 
the beginning of the observing time, the first point in the 
X-ray light curve lies at ~ 2.5 a from the model. The X-ray 
behavior might also be consistent with the optical/UV light 
curve (Arcavi et al. 2011). A bump in the X-ray light curve, 
even if not statistically significant, can also be observed at 
~ 13 d after the SN explosion. This is barely consistent with 
the peak of the Swift/UVOT light curve in the U band, when 
the SN peak is clearly visible (Arcavi et al. 2011). 

We extracted two set of spectra. The first comprise the 
entire dataset. Assuming the same absorption pattern de- 
rived from XMM-Newton observations we can fit the 414 
counts in the 0.3-10 keV energy band (wit only 71% due 
to the source) with a power law model. The best fit pho- 
ton index is F = 1.48 ± 0.18. The reduced x 2 cd = 1-5 with 
18 d.o.f. (7% null hypothesis probability). The average flux 
0.3-10 keV unabsorbed flux is 1.2 x 10" 13 erg cm" 2 s" 1 
(giving a luminosity of 1.0 x 10 39 erg s _1 ). A double MEKAL 
model improves only marginally the fit. This is expected be- 
cause spectral evolution was observed, at least, across the 
two XMM-Newton observations. A spectrum was also ex- 
tracted taking data from the first five days, summing 30 ks 
of data and collecting 134 photons (5% background). In this 
case the spectrum is harder, with a power law photon index 
of T = 1.1±0.3 (reduced Xred = L2 for 11 d -°- f -. 27% n.h.p.) 
The 0.3-10 keV unabsorbed flux is three times more than 
the average value. 



4 DISCUSSION 

X-ray absorption holds the potential of tracing the pres- 
ence of metals through their impact on the X-ray spectrum. 
XMM-Newton observations clearly indicate a non-negligible 
absorbing column density in addition to the Galactic con- 
tribution. This likely arise within M51. Assuming a solar 
distribution of elements and a solar metallicity the column 
density is logAu = 20.7l o ;^ (90% confidence level). This 
value can be compared to optical estimates. Arcavi et al. 
(2011) estimated E(B - V) £ 0.05. This can be converted 
to log ~ 20.5, following Predehl & Schmitt (1995). Data 
from the The H I Nearby Galaxy Survey (THINGS; Walter 
et al. 2008) yields at the position of SN 2011dh log N(H I) = 
20.6 ± 0.1. A more thorough analysis was carried out by 
Ritchey & Wallerstein (2012). From echelle observation of 
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[htb] 

Table 2. XMM-Newton spectral fits. 



Model 


Column density 
(10 20 cm 2 ) 


1 st Obs. 
T/feT (keV) 


2 nd Obs. 
T/kT (keV) 


x? ed / dof 

(n.h.p.) 


Power- law 
Black body 
Bremss. 
Mekal 


< 0.9 

o 7+ 2. 3 
z - '-2.1 
2 3+ 21 


I. 6 ±0.2 

u ' oa -0.04 

II. 9± 25 7 5 


2.0 ±0.3 

0.34+To2 
2.8l?;° 2 

3 0+ 24 


1.4/58 (0.02) 
3.4/58( 10- 16 ) 
1.5/58( 0.01) 
1.4/58( 0.02) 


Mekal+pow 


67.2+ 27 6 5 


0.22 ± 0.01 
1.0 ±0.2 


29+ - 45 
U ' ZM -0.15 

1.3 ±0.4 


1.3/54 (0.07) 


Mekal+Mekal 


4 7+5-5 


> 10.7 


> 2.4 
59+ 26 


1.2/54 (0.15) 



Errors have been calculated for a A% 2 = 2.71, corresponding to a 90% confidence level for one parameter of interest. 
The column density reported in the Table is in addition to the Galactic column density {Njj = 1.8 X 10 20 cm -2 ) and it likely is within 

the host galaxy. 

Leaving free the metal abundance of the MEKAL model, does not improve the fit, obtaining only an upper limit on the metallicity of 

Z < 3.9 Zq. 




Time (d since explosion) 



Figure 3. Swift XRT light curve starting from the best estimate 
of the SN explosion time. The light curve has been obtained using 
the Leicester University user_object tool, by adaptively summing 
up data to reach a 6 a detection. The two diamonds indicate the 
times of the two XMM-Newton observations. 

the SN, they derived Na I column densities of the multi- 
ple components and estimated log N(H I) ~ 19.8 — 20.2. 
The value derived from X-ray observations is slightly larger 
than this but still consistent. This would indicate that the 
surrounding medium has a metallicity similar to the solar 
value. 

One question is if this instrinsic absorbing column den- 
sity is intrinsic to the source or due to some intervening 
matter in M51 along the line of sight but far away from 
the SN site. M51 is a nearby galaxy and it is well studied 
in the X-ray band. M51 hosts a number of ultra-luminous 
X-ray sources (ULXs) well studied and characterized (De- 
wangan et al. 2005; Yoshida et al. 2010). The two closest 
ULX sources to SN 201 ldh are source 82 (~ 1 arcmin to the 
north) and source 63 (~ 1.5 arcmin to the north-west). Both 
sources are variable. The column density to source 82 is at 
a level of log An ~ 21 (Yoshida et al. 2010). We confirm 



this value with the analysis of our XMM-Newton data (as- 
suming a power law spectrum). Source 63 is fainter and the 
intrinsic column density (assuming a power law plus MEKAL 
model) is consistent with zero (Dewangan et al. 2005). In 
our data source 63 is a factor of ~ 3 fainter and we can- 
not confirm these findings. We also note that the Type Ic 
SN 19941 occurred 2.3 arcmin away from SN 2011dh. No 
direct evidence for the presence of circumstellar matter was 
found in X-ray data (namely too few photons to assess an 
increase over the Galactic column density value), but the 
X-ray emission was interpreted as due to the interaction of 
the blast wave with its surrounding circumstellar material 
(Immler, Wilson & Terashima 2002) . From this comparison 
we can assess that we have, at least, found strong indications 
that the intrinsic column density observed in SN 201 ldh is 
due to circumstellar material, as expected for a supergiant 
progenitor (Bersten et al. 2012). 

X-ray emission from SNe lib has been observed in a 
very limited number of objects: SN 1993J (Chandra et al. 
2009), SN 2001gd (Perez- Torres et al. 2005), SN 2001ig 
(Schlegel & Ryder 2002), SN 2003bg (Soderberg et al. 2006), 
SN 2008ax (Roming et al. 2009), and SN 2011dh. SN 1993J is 
probably the best studied object. Using Rosat data, Zimmer- 
mann & Aschenbach (2003) found that the intrinsic column 
densities is much higher than the Galactic column density. A 
deeper look came from ASCA data where Uno et al. (2002), 
fitting a two-component thermal-plasma model, found that 
the low-temperature component has much higher column 
density than the high-temperature component. These high 
column densities in excess to the Galactic absorption were 
attributed to the absorption by an additional cool shell. For 
all the other X-ray detected SN lib we have too few counts 
to assess if an intrinsic column density excess is present. 
Thanks to the larger XMM-Newton throughput, we were 
able to assess that the X-ray spectrum consists of one soft 
and one hard emission components (modeled with the MEKAL 
model), likely related to a reverse shock as in SN 1993 J 
(Chandra et al. 2009; Nymark et al. 2006). Our best fit two 
component thermal plasma model is then fully consistent 
with shock model physics from SN explosions, i.e., a hard 
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thermal plasma component (at a few keV) that arises from 
the fast forward shock into the circumstellar material, and 
a soft thermal plasma component (below 1 keV) from the 
emergence of a reverse shock at lower speeds into the ejecta. 
Since the reverse shock is located behind the forward shock, 
it is expected to have a higher absorbing column density 
than the forward shock. Our XMM-Newton observations 
clearly show a non-negligible absorbing column density in 
addition to the Galactic contribution. 

Finally, a reanalysis of the Swift XRT light curve might 
indicate a flux excess at the beginning of the X-ray obser- 
vations. Despite the fact that the statistical significance of 
this increase is low and the physical interpretation is chal- 
lenging, this excess is likely due to the adiabatic cooling of 
stellar's photosphere a few days after the shock breakout. 
The time of the X-ray excess is identical to the time of the 
optical excess described by Arcavi et al. (2011) and inter- 
preted by them as the tail of the shock breakout emission. 
This might hint for a non-compact progenitor, i.e. yellow 
supergiant progenitor. 
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